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Abstract 



'i> We present an interferometric technique for the reconstruction of ultra-wide band impulsive sig- 

^n nals from point sources. This highly sensitive method was developed for the search for ultra- 

^^ high energy neutrinos with the ANITA experiment but is fully generalizable to any antenna array 

ly-v detecting radio impulsive events. Applications of the interferometric method include event re- 

construction, thermal noise and anthropogenic background rejection, and solar imaging for cal- 
ibrations. We illustrate this technique with applications from the analysis of the ANITA-I and 
ANITA-II data in the 200-1200 MHz band. 
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H 1. Introduction 



In the last decades there has been an increased interest in using radio frequency (RF) instru- 
mentation for the detection of ultra-high energy (UHE) > 10'^ eV neutrinos and cosmic rays. 
Russian- Armenian physicist Gyurgen Askaryan |[T] predicted that high energy particle showers 
produced in dense dielectric media would result in impulsive coherent Cherenkov radiation at 
radio frequencies. The emission was experimentally confirmed for the first time in 2001 using 
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showers induced by high energy photons in sihca sand f2l. The results are consistent with mod- 
ern particle shower and radio emission simulations [3J ^nd this effect has since been observed in 
salt H and ice 0. Several experiments exploit this technique in the search for UHE neutrinos 
using antennas buried in ice fS" 7|, radio telescopes pointed at the Moon |l8l|9l, or balloon-borne 
antenna arrays orbiting the Antarctic continent |10|. 

Cosmic ray extensive air showers (EAS) produce a radio impulse due to the transverse current 
produced by the separation of electrons and positrons resulting from interaction with the Earth's 
magnetic field I.11..12J . This geo-synchrotron emission was first observed in a ground array in the 
1960's ifTSl . Since then, there have been many observations lfT4l - ll25l . Recently, the Antarctic 
Impulsive Transient Antenna (ANITA), a balloon-borne antenna array that synoptically scans 
the Antarctic continent in the 200-1200 MHz range, observed geo-synchrotron emission in the 
ultra-high energy range for the first time ll26l . 

The growth of this field demands improved analysis techniques. In particular, it is expected 
that the first neutrino observations will be from signals that are close to the detector threshold set 
by thermal noise. This requires analysis techniques that are highly sensitive and can efliciently 
discriminate between a weak impulse and a thermal fluctuation. 

Interferometric methods have been widely and successfully used in radio astronomy. Radio 
telescopes are able to map weak sources in the sky using the correlations between signals in an 
antenna array. Distant sources are imaged via the relation between phase delay and source direc- 
tion |27|. Additional point-spread deconvolution methods are applied to reveal high resolution 
brightness maps of the radio source. The most precise astrometric measurements (200 /i-arcsec 
resolution) are obtained from the correlations of a single pair of antennas with 8,000 km separa- 
tion 1,28 J . 

The signals produced by UHE particle showers are rather diff'erent from those detected by radio 
interferometric telescopes. UHE particle shower emissions are impulsive transient events (on 
microsecond to nanosecond time scales) while radio astronomical sources are better described 
as brightness distributions that can be imaged with long exposures. Despite the differences in 
the nature of these signals, the fundamental ideas of radio interferometry can be applied to the 
detection and analysis of UHE particle shower impulsive transient events. 

An interferometric approach for impulsive signals was developed for ANITA and successfully 
applied to the data analysis of both flights ||26l|29]| . The ANITA antenna array is designed to 
observe impulsive radio emission in the frequency range 200-1200 MHz from UHE neutrinos 
interacting in the Antarctic ice. Each antenna is a dual-polarized quad-ridged horn with a gain 
of ~10 dBi and a half-power beam-width of 45°. The ANITA array is cylindrically symmetric 
with neighboring antennas that have a center-to-center distance of 1 meter. Typical antenna 
separations used in the direction reconstruction are 5 meters for ANITA-I IfTOl and 7 meters for 
ANITA-II 1 29 1 producing a pointing error below 1 ° . 

In the interferometric imaging procedure developed for ANITA, each pair of antennas pro- 
duces a fringe across the sky at an angle corresponding to the baseline directiorTl The fringes 
are then summed together resulting in an image that peaks at the source location. Such an image 
is named a "dirty map" in radio astronomical usage, and reduction of sidelobes is possible with 
further image processing. However, we have found that for this type of "pulse-phase interfer- 
ometry" HIOL these maps are adequate since the sources of interest to ANITA are unresolved. It 
is also worth mentioning that ANITA does not apply this mapping to identify and characterize 



A baseline is the vector defined by an antenna pair. 



source structure but rather for the identification of coherent point source impulses. We also rely 
on the point-like characteristics of our data and delay closure to calibrate antenna positions and 
cable delays |10|. 

Although there have been other impulse beam-forming results in the past, particularly the one 
developed for imaging the radio flashes from ultra-high energy cosmic rays (UHECRs) for the 
LOPES antenna array [i30l, there are some important differences with the variant developed in 
this paper ANITA is a self-triggering array and does not have muon counter data for identi- 
fying UHECR signals. The interferometric techniques developed for ANITA are applied as a 
stand-alone technique for identifying plane wave impulses in the data and for refining the pre- 
cision of the directional reconstruction while providing improved rejection of thermal noise and 
anthropogenic backgrounds. 

This paper presents an interferometric method applicable to broadband, impulsive radio sig- 
nals for antenna arrays with large fields of view and digital waveform recording capabilities. The 
technique is illustrated via its application to the ANITA analysis. In Section 2 the mathemati- 
cal foundations of the interferometric image production applied to radio impulses are covered. 
Section 3 describes the application of interferometric images to point source impulse reconstruc- 
tion along with rejection of thermal noise and anthropogenic backgrounds. Section 4 presents 
an application of the interferometric method to identify and characterize sources that are below 
detection threshold but continually present. We demonstrate the technique with observations of 
the Sun and images of RF activity on the Antarctic continent. In Section 5 we conclude this 
paper and mention some future applications of the interferometric technique developed here. 



2. Interferometric Equation 

In this section we formulate the interferometric approach used for radio impulses. At its core, 
interferometry is based on combining multiple measurements of the same signal. We discuss the 
relation between a recorded voltage and an electric field followed by the relation between the 
Adding Interferometer and a Cross-Correlation Interferometer in the context of radio impulse 
detection. We motivate the approach used for the analysis of ANITA data. It is worth noting, 
that unlike typical interferometric arrays, the ANITA antennas are not all pointing in the same 
direction as their boresight direction varies with payload azimuth to provide a full 360° field of 
view coverage ifTOl . In general, this would require that the system impulse response be decon- 
volved prior to beam-forming. However, we argue that for the ANITA horn antennas this is not 
necessary. The tools described below can just as well be applied with prior deconvolution of the 
antenna response but not requiring this step makes their application practical. 

2.1. Relation between the electric field and receiver voltage 

An incident electric field couples to the antenna inducing surface currents that produce voltage 
differences in a transmission line that can be stored in a recording apparatus. Since ANITA is 
primarily concerned with ultra-wide band impulses it is natural to approach the problem in the 
time-domain. A detailed time-domain treatment of the relation between an electric field E(f) 
and the voltage v(t) recorded in an experiment can be found in Ii3 1 J and references therein. The 
relation is captured via 

v(f,f)= .&h.v.(f)*E)(0 (1) 
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where Z^ is the load impedance, Zq is the impedance of free space, and f is the direction of 
incidence of the radiation. The effective height vector h5,.j(f , t) is the time-domain representation 
of the antenna receiver system complex impedance, and is equivalent to the antenna receiver 
system response to a delta-like pulse 1 3 1 1 . The • operator is a vector convolution defined by 
(f * g)(f) - J dst(s)-g(t - s). We have made hjvj(f , t) explicitly a function of f since, in general, 
the effective height depends on the angle of incidence of the electric field with respect to the 
antenna. The measured voltage v recorded with a digitizer is, strictly speaking, a function of 
time only; the added dependence of v on f captures the fact that the frequency contents and 
group delay of the effective height depends on the direction of incidence of the radiation relative 
to the antenna. 

Figure [T] shows the time-domain effective height for various incidence angles for an ANITA 
quad-ridged horn antenna. It is important to note that Figure [T] shows the effective height of 
the antenna alone and does not include the system response of the full signal chain with cables, 
filters, and amplifiers. The low frequency dispersion (corresponding to the frequency range 
200-300 MHz) extends for about 10 ns on the tail end of the waveform while the high frequency 
(300-1200 MHz) portion the signal is contained within the first few nanoseconds. Note, however, 
that the effective height dispersion is nearly identical for signals up to 45° away from boresight. 
This means that, for the purposes of correlating signals, the antenna response does not need to 
be corrected within this angular range. The only significant difference in the effective height 
function, for the various observation angles within 45° away from boresight, is the attenuation of 
high frequencies. This feature is equivalent to the antenna beam pattern being narrower for high 
frequencies and wider for lower frequencies. 
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Figure 1 : ANITA quad-ridge hom antenna eflfective lieight for vertical polarization. The effective height is plotted for 
various elevation angles. The effective height near boresight has the strongest response at high frequencies showing the 
highest peak. As the incident electric field is moved away from boresight, the effective height loses sensitivity to the 
highest frequencies but retains sensitivity for the lowest frequencies. The ringing in all of these waveforms is dominated 
by the low frequencies (< 300 MHz) and is dispersed over tens of nanoseconds. The similarity of the antenna response 
over a wide range of angles around boresight allows for interferometric reconstruction of source positions without the 
need of deconvolving the signals with the directionally dependent efl'ective height function. 



2.2. Mapping the incident signal direction using receiver voltages 

The similarity between waveforms due to the same signal in the antenna array is at the 
foundation of interferometry. Radio telescope interferometers typically use the pairwise cross- 
correlation between signals recorded at each antenna where the waveforms are delayed according 
to a given direction and multiplied together. Another variant is to use the Adding Interferometer 
where the waveforms of each antenna are delayed, summed together, and the square of the sum 
is integrated. Summing the waveforms, delayed according to the direction of incidence, averages 
down the noise while coherently adding the signal, providing an improved signal to noise raticPl 
(SNR). This summed waveform has the clear advantage of exposing a weak signal measured by a 
number of antennas. In this section we derive various interferometric quantities from this starting 
point. Although the quantities associated with the Adding and Cross-Correlation Interferometers 
are well known in the literature |27|, we re-derive them here in the context of radio impulses to 
motivate the analysis techniques used in the next section. 

Let us describe the data of an array of A^^i antennas by a set of voltages v,(r) for each antenna 
;. For a plane wave incident from direction f the delay T,(f) at each antenna is given by 

T,(f) = - (D - R; ■ f) (2) 

c 

where D is the distance between the source and the antenna array, R, is the antenna position, as 
shown in Figure [2] and c is the speed of light. The phase aligned voltage waveforms are to be 
summed together to give the coherently summed waveform 

Vs(f,f) = 2v,(f + T,(f)). (3) 

1=1 

If the voltage waveforms v, are only due to uncorrected noise between the antennas, with the 
same root-mean-square (RMS) voltage v„„.5, then V^. has an RMS increased by a factor of the 
VA^, regardless of the delays between them. If the voltage waveforms all contain the same plane 
wave signal impulse, then Vy, will be equal to Na times v,- when the delays correspond to the 
direction of incidence of the signal. Thus, on average, a set of waveforms with the same noise 
RMS and the same signal will result in V^ with an amplitude enhanced by a factor VA^ over 
each individual v,. 

Figure [3] shows a model of the ANITA horn antenna array with ten signals highlighted. The 
event shown is sent from a ground-based calibration pulser used for testing pointing reconstruc- 
tion techniques. The top panel of Figure HI shows the phase (or delay) aligned waveforms using 
the known direction of incidence, which are summed to give the coherently summed waveform 
Vj, shown on the bottom panel of Figure p] 

One could formulate an analysis based on finding the peak of Vy, as a function of f , but in 
reality signals are not perfect delta functions and display some dispersion due to the antenna, 
the signal chain, or the neutrino induced shower itself |32l. For this reason, one can obtain a 
higher SNR measurement using the power of Vy, integrated over a time window T relevant to the 
impulses of interest. 



^We define the signal to noise ratio of an impulse by the half maximum peak to peak voltage difference divided by 
the root mean square of the noise. 




Figure 2: Geometrical delay for a pair of antennas. The delayed combination of waveforms is the foundation of interfer- 
ometric mapping. The schematic above shows how the delay of two signals is dependent on the geometry. Vectors R; 
and Ry denote the position vectors of the two antennas. For an electric field coming from the incident direction f from 
a distance D the geometrical delay Ar between the two antennas in the far-field is given by cAr = f • (R, - Rj) where 
c is the speed of light and R, - Ry is the baseline vector. Note that the distance of the source D is not present in the 
far-field approximation of the delay. ANITA events come from >100 km away for meter scale wavelengths justifying 
this approximation. 

The time averaged power of the summed receiver vohages is 

i^s(f)=^7 r dtVliur) (4) 

^L i Jo 

where T is the total time of the integration and Zl is the impedance of the system. The quantity 
f 2(f) is also known as the Adding Interferometer. 

FigurelSlshows the power map Pz{f) for a flight calibration impulse from the ANITA-I flight. 
A set of ten antennas centered around each 0-sector] shown in Figurepl is used for the summed 
waveform as a function of direction. The power map shows the most likely direction of incidence 
as a large peak with sidelobes representative of the system's point spread function determined 
by the geometrical arrangement of the antennas and the interference pattern of the waveforms. 
At angles away from the direction of incidence of the impulse, the image shows the typical 
random pattern produced by thermal noise. Note that this approach has ignored the differences 
in effective height between antennas pointing in diff'erent directions. This is because the ANITA 
horn antenna response does not significantly vary in phase at different incident angleq^ It also 
allows for faster computation which is advantageous when dealing with large data sets. 

In the context of mapping brightness distributions, the Adding Interferometer has the poten- 
tially undesirable feature of including the noise from each individual waveform. It also has 



'ANITA is divided into 16 (^-sectors, each consisting of a pair of antennas on top and on bottom (see left side of 
Figure[3l. ANITA-II has a third tier of antennas on every other ^-sector 

''Finite Difference Time Domain simulations of the impulse response of the ANITA horns found 45 ps delays between 
signals incident on boresight and at 22.5° away. At 2.6 Gsa/s the digitization time bin width is 384 ps. 
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Figure 3: A schematic of the ANITA horn antenna array. Each antenna is a dual polarized quad-ridged horn. The ten 
antennas closest to the direction of the incident impulse are highlighted with the recorded waveforms shown on the right. 
The signal is from a ground-based calibration pulser and only the vertically polarized channels are shown. The additional 
ringing in these impulses is due to a combination of the antenna response (see FigurefTJ in addition to the filters in the 
signal chain and the ringing of the transmitted impulse itself. Note that the signals are very similar within the five phi- 
sectors highlighted. The main difference is in the reduced high frequency response of the impulses detected away from 
boresight which is due to the directional dependence of the antenna effective height shown in FigurefT] 



undesirable effects when the gain and noise figures in each channel are not matched. Although 
P^{r) can be a useful quantity for the analysis of radio impulses, given that the proper calibra- 
tions have been made, we can also remove some of its potentially inconvenient qualities with 
some mathematical manipulations. 

The time averaged power of the summed receiver voltages in Equation!?] can be expanded to 
give 



Pz(r) 



J J Na Na ^t 



Jo 



(t + Ti(r))vj(t + Tj(r)) 



(5) 



where t, and tj are the delays with respect to the origin of a given coordinate system in Equa- 
tion|2] The terms on the right hand side of the equation are proportional to the cross-correlations 
between antenna voltages defined as 



V; ® Vj(f) = I dt Vi(t)Vj(t - ATij) 

Jo 



(6) 



where At, 



ij - '' 



Note that the term D, in Equation l2] denoting the distance from the source 



to the array vanishes. Substituting Equation [6] into EquationjS] reduces to 



PAf) = J]Pi 



1 1 



Na 

vv 



^ 1=1 Itl 



'i ® Vj(r) 



where 



1 1 r 

Pi^—- dtvj(t) 

^L I Jo 



(7) 



(8) 



is the average power of each individual waveform, which does not depend on f given that t,-,- - 0. 
The cross-correlation (also known as cross-power) term contains the directional information of 
the electric fields incident on the antenna array. 
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Figure 4: Top: phase aligned voltage waveforms from Figurelslaccording tlie to tlie known direction of incidence. Given 
the positions of the ground-based calibration pulser and the payload, the direction of incidence in payload coordinates 
can be determined with a high level of accuracy. This information is used to translate the direction of incidence to the 
geometrical delays expected on each antenna. Bottom: the coherent sum Vj;(f, f ) of all ten waveforms above. 



If we only keep the terms that depend on the direction of incidence f in Equation 17] we obtain 
cross-correlation map 



Z 



"/I 



(9) 



where the restriction j < i counts each baseline once. However, the quantity X(x) retains unde 
sirable features if the gains and noise figures of each channel / are not matched. 

Another approach is to normalize the cross-correlation by the power of the waveform accord 
ing to 

V; ® Vj(f) 



Cum 



ylj,^dt\Vi(t)\2yljJdt\vj(t)\2 



(10) 



which is known as the the cross-correlation coefficient or coherence function fSSl. This value is 
bounded between a maximum value of -i-l and a minimum value of -1 and quantifies the similar- 
ity between waveforms / and j. If the waveforms are identical the cross-correlation coefficient is 
-1-1. If they are identical with a 180° phase difference then it is equal to -1. The more dissimilar 
the waveforms, the closer the cross-correlation coefficient is to zero. 
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Figure 5: Top: Map of the time averaged power as a function of incident direction Pj:(r) for the signals shown in 
Figurel4] The peak of the image corresponds to the direction of incidence of the radiation. However, rather than mapping 
the source structure of an object, this image identifies the location of a point source impulse. The sidelobes are due to the 
geometrical arrangement of the antennas. Bottom: azimuth and elevation slices along the peak of the power map. The 
azimuthal slice (left) shows the summed waveform power both in the region of the peak and on the opposite side of the 
payload where only thermal noise is present. The power levels are consistent with system noise temperature of 320° K. 
The elevation slice (right) has an offset due to all antennas along a phi-sector being pointed in the same direction. The 
minimum level of 3 pW rcT^ is determined by the summed power of each antenna while the oscillations are due to the 
interference terms. 



The power sum can be written in terms of the cross-correlation coefficients as 



(11) 



i*J 



M this sense the cross-correlation coefficients naturally quantify the interference terms in the co- 
herent power sum. The C,j(f) terms are not sensitive to the overall amplitude scale of the wave- 
forms, which can vary due to thermal fluctuations. The Cy(f) terms contain all the directional 
information provided by Pz{r) of the Adding Interferometer or X{f) of the Cross-Correlation 
Interferometer 

For the ANITA analysis it was found that the cross-correlation coefficient was the best means 
of reconstructing the direction of a signal |10l [34l [26] |29l. The attractive features are that it 
normalizes out overall amplitude fluctuations along with mismatched gains and noise figures. In 
addition, its statistical behavior is not strongly affected by the moderate use of notch filters. 

Another way to produce interferometric images is to project of the cross-correlation coeffi- 
cients Cij{r) of antenna pairs onto the incident angle space. We define the coherence map of a 
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set of waveforms as 



1=1 j<i 



(12) 



where the restriction j < i is put in place so as not to count the contribution of each baseline twice 
and Nb is the number of baselines formed by a set of A^a antennas given by Nb - Na{Na - l)/2. 

Figure l6] shows the projection of a cross-correlation coefficient from signals incident on a pair 
of antennas for several orientations. The time-domain fringe pattern is projected to the incident 
direction space, in payload elevation and azimuth coordinates, and their direction is perpen- 
dicular to the baseline vector orientation. The azimuthal resolution is dominated by horizontal 
baselines while elevation angle resolution is dominated by the vertical baselines with contri- 
butions from the diagonal baselines. The fringe width in the incident direction plot is inversely 
proportional to the separation of the antennas and depends on the frequency content of the signal. 

An example of the image formed by the coherence map M(f ) is shown in FigurelT] The image. 
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Figure 6: An example of the formation of an interferometric image. From top to bottom tlie cross-correlation coefficient 
of an antenna pair is displayed. On the left, the antenna pair used for the cross correlation is highlighted. In the middle, 
the cross-con'elation coefficient of the waveforms corresponding to the antenna pair on the left is shown as a function 
of time. The signals in Figure[3]were used and they display a strong cross-correlation coefficient. The section of the 
waveform highlighted in red corresponds to delays whose geometry is constrained by the field of view of the antenna 
pair. On the panels to the right, the same cross-con'elation coefficient is plotted as a function of incident direction in 
payload elevation and azimuth coordinates which is related to the delay via Equation 3. The true incident direction of 
the radiation lies at an elevation of -8 degrees and an azimuth of 259 degrees. Note that the directional projections of the 
cross-correlations all overlap at this point. 
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formed by the superposition of fringes oriented in various directions, peaks in the direction of 
incidence of the impulse. Note that although the individual cross-correlation coefficients C,y(f ) 
have sidelobes comparable to the true direction of incidence, their superposition produces a sharp 
peak, which greatly reduces the possibility of mis-reconstruction and increases the ability to 
reconstruct the direction of noisy signals. 

In general, a windowing strategy is necessary to produce a coherence map M(f) and the choice 
depends on the properties of the antennas as well as the geometric configuration of the array. In 
the case of ANITA, the coherence map M(f ), as shown in Figure IT] and subsequent figures, is 
based on the symmetry of the array as well as the antenna beam width. The global image is 
composed of the stitching of 16 images, each with an azimuthal field of view of 22.5°, centered 
around each (^-sector of the pay load. All the antennas from four adjacent 0-sectors ( 5 0-sectors 
in total) are used in the formation of the image. 

Figure Islprovides examples of the coherence map M(f ), coherent waveform sum V^{t, f) in the 
direction of the peak of M(f ), and its power spectrum for signals from the ANITA-I flight. The 
cosmic ray signal, shown on the left, is strong with a clearly defined peak on the coherence map. 
Although the cosmic ray signal is highly impulsive fSSl, the coherent waveform sum displays 
various oscillations due to the low frequency contents of the signal, where the ANITA impulse 
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Figure 7: Top: A global coherence map M(f) u.sing all antenna pairs sharing the same field of view from Figure|6] The 
full field of view image is a composite of coherence maps M(f) with an azimuthal span of 22.5° centered around each 
(/i-sector (see text for details). The azimuthal slice along the peak (bottom left) has very small values of M(r) away from 
the peak since all the antennas forming the image are seeing thermal noise. The elevation slice along the peak (bottom 
right), however, shows a larger degree of cross-correlation away from the peak. This is due to the fact that all antennas 
are pointed in the same direction. The formation of this image is the "dirty map" of radio astronomy. 
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response is dispersive (see Section 2.1). In the center of FigurelS] a 450 MHz carrier-wave (CW) 
signal displays a coherence map with multiple lobes of comparable amplitude. The 1150 MHz 
signal, shown on the right, displays many more lobes of smaller angular size on the coherence 
map. 

In this section we have derived the Adding Interferometer equation for Pzir) and Cross- 
Correlation Interferometer equations for X{r) and M(f ) from the coherent waveform sum V^it, f ). 
Although we primarily use M(f ) for the ANITA analysis, this is not to say that X{r), P^ir), and 
V^(t, f ) do not have advantages for certain applications where it makes sense to use them. As we 
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Figure 8: Examples of coherence maps M(f ) for different types of signals. From left to right a cosmic ray impulse, a 
450 MHz carrier wave (CW), and a 1150 MHz CW signal are shown. In the middle row from the top. the coherently 
summed waveforms corresponding to the peak of the image are shown. The cosmic ray impulse shows the dispersion 
characteristic of the ANITA signal chain at low frequencies. The carrier waves display a strong sinusoidal behavior. On 
the bottom plots the power spectral densities are shown for the direction of the peak of the coherence map. The cosmic 
ray impulse is broadband with a falling spectrum while the carrier wave signals are strongly peaked at a single frequency. 
The thermal noise power spectral density is shown in red. 
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will show in the next section, the ANITA analysis relies heavily on the combination of M(r) with 
the peak of V-^(t, f). With the relations between the interferometric quantities presented in this 
section being well understood, an analysis can be tailored to use each quantity, or combinations 
of them, as needed for the desired application. 



3. Reconstruction of Impulsive Point Sources 

The coherence map M(f ) is a powerful tool not only because it allows for the identification of 
a point source direction but also because it provides an efficient means to reject thermal noise and 
weak sources of interference. The coherence map can be combined with the coherent waveform 
sum Vz{t, f) and other related quantities from the previous section to provide for more sensitive 
analysis tools. In this section we discuss point source reconstruction and background rejection. 
We illustrate the application of the tools developed in this paper with examples from the ANITA 
analysis. 

3.1. Event Reconstruction 

The primary means of reconstructing a signal is to identify the peak of the coherence map 
M(f). The coherent waveform sum Vz(t, f ) provides a means of inspection for the quality of 
the reconstruction. One can compare the peak of the coherence map to the peak of yj;(f, f) or 
P^ir) as an additional check of the correctness and fidelity of the reconstruction. The individual 
waveforms v,(f) should appear aligned, as was shown in Figure l3] 

For an interferometric image, such as M(f), the resolution of a point source is proportional to 
A/R, as prescribed by the Rayleigh criterion, where A is the wavelength of the radiation and R 
is the separation between antennas. It is important to distinguish between the resolution of the 
image, which is the width of the peak, and the reconstruction error, which is the statistical scatter 
on the location of the peak after repeated measurements. In the case of wideband radio impulses 
the errors can be significantly better than the resolution estimate provides. The factors that affect 
the error of the location of the peak, besides the ratio A/R, are the SNR, the bandwidth of the 
impulse, and the number of antennas observing the signal. 

For uncorrected noise, the error on the source direction is inversely proportional to SNR (see 
Appendix A). The fact that a wideband impulse does not reside in a single frequency but rather 
over many bands requires some care in relating the resolution (proportional to A/R) to the signal 
to noise ratio. The relation between the SNR of an impulse and the signal-to-noise ratio of its 
spectral components depends on the spectral shape of the impulse and the noise background. 
For a digitizer with sampling frequency /, the frequency resolution is A/ = fs/N^ where A^^ is 
the number of samples recorded. For a signal with bandwidth B, the number of independent 
frequency measurements made on the signal is Nf = B/Af. For a non-dispersive signal with a 
constant spectral signal to noise ratio (where the signal Fourier amplitude is a constant multiple 
of the average thermal noise Fourier amplitude for a given frequency range), the time-domain 
impulse SNR is related to the signal-to-noise ratio of the individual spectral (snr) amplitudes by 
SNR = snr ^JBTKf, where B is the bandwidth of the impulse. In this illustrative case, the error on 
the direction of a wideband impulse is proportional to the A/R, where A is the central wavelength 
and inversely proportional to snr yJBjXf. See Appendix A for a more rigorous derivation. 

The improvement in directional reconstruction errors can be estimated by counting the number 
of independent measurements contributing to the result. In the paragraph above we discussed the 
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contribution of the number of bands B/Af. For an array with A^^ antennas, there are A^^ - 1 in- 
dependent baselineaj The number of independent measurements contributing to the directional 
reconstruction of the signal is the number of independent baselines times the number of indepen- 
dent frequency measurements. For the ANITA antenna geometry, the directional reconstruction 
error improves approximately by a factor of V(-^a - 1) over the image resolution estimate pro- 
portional to A/R. See Appendix A for a more rigorous estimation of the point source impulsive 
reconstruction errors. 

For the ANITA analysis, the coherence map M(f) is calculated with 1° x 1° pixels. The peak 
widths of impulsive events for ANITA are typically ~ 3° for elevation and ~ 10° for azimuth 
(see Figure ITli depending on the frequency contents of the impulse (see Figure IHll. The peak 
widths are consistent with the A/R resolution estimate for the frequencies and basehne lengths 
involved in making the image. When testing the angular reconstruction errors for point sources 
with a calibration pulser, the statistical errors are measured to be 0.26° in elevation and 0.56° in 
azimuth [29.1 . For ANITA, the sampling frequency /, - 2.6 Gsa/s and the recording window is 
256 samples giving A/ ~ 10 MHz. The calibration impulse has a bandwidth B ~ 300 MHz. 
Each pixel in the coherence map image uses 10 antennas. The angular error therefore improves 
by ^j{N~A^^yyB]~KJ ~ 16 over the resolution, resulting in an estimated errors of 0. 18° in elevation 
and 0.5° in azimuth. The elevation and azimuth errors are consistent with the measured errors of 
0.26° and 0.56°, respectively. 

Both the ANITA-I and ANITA-II data show a clear dependence on angular error with SNR. 
Figure [9] shows the dependence of angular error on single antenna SNR for the ANITA-II point- 
ing calibration. The quality of the angular reconstruction error can vary by a factor of 2 between 
the weakest and strongest signals. The comparison between the detected SNR, plotted in Fig- 
ure l9] and the theoretical SNR used in the discussion above is not straightforward due to various 
complications such as impulse dispersion (the signal is not a perfect delta-function impulse), the 
impulse signal and noise spectral snr varies with frequency, and there are additional error con- 
tributions due to clock synchronization and other cahbrations. However, the discussion above 
does provide a reasonably accurate estimate of the pointing error improvement over the image 
resolution. 

3.2. Thermal Noise Rejection 

The ANITA data consist of >99% thermal noise events with potentially a few neutrino events 
expected near the thermal noise threshold. The ANITA analysis is searching for a small signal 
sample in a large thermal noise background and therefore requires a highly efficient thermal noise 
filter. This section will discuss how the interferometric quantities developed in this paper were 
used to successfully obtain a highly sensitive means of rejecting thermal noise. 

The combined use of the coherence map Mir) and the coherently summed waveform Vj,{t, f ) 
provide a complementary handle on thermal noise rejection. Figure[TO]shows the peak coherence 
map value M{r„,ax), where r„ax is the direction of the peak, plotted against the SNR of VY.{t, r„ax)- 
In the case of thermal noise, these quantities become anti-correlated when either of them have 
a high value. If the total level of coherence fluctuates upward, the peak value of the waveform 
sum tends to be low. Conversely, if the noise displays a large peak, the correlation, which uses 
all points in the waveforms, tends to be low. This is an advantageous discriminator of impulsive 
signals which have a high peak coherence map value and a high summed waveform SNR. 



The total number of baselines Na (Na - 1 )/2 can be represented as linear combinations of a subset of Wa - 1 linearly 
independent baseline vectors. 
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For this purpose we apply a Fisher Hnear discriminant, in the form of L = y+mx, that combines 
the values of the peak coherence map y and coherently summed waveform SNR x. The slope 
m of linear discriminant is chosen so that it is tangential to the contours shown in the thermal 



noise simulations in Figure 10 while L shifts the overall level up and down. The distribution of 
the linear discriminant value L, shown in Figure [TT| has an exponential fall-off at the tail. The 
trend has been extrapolated to set a cut consistent with 0.01 thermal events passing for all of 
the ANITA-I data set. The value of this discriminant is shown for the all ANITA-I events in the 



vertically polarized channel. The black dots in Figure 10 show the distribution of 1000 simulated 



triggered neutrino signals ifTOl showing that the linear discriminant is a highly efficient filter of 
thermal noise. 



3.3. Radio Frequency Interference 

Prior to the pointing analysis, the ANITA data goes through an adaptive notch filter to iden- 
tify and remove CW peaks from signal spectra of the data. The most common source of mis- 
reconstruction is the presence of a weak CW signal that has survived the pre-filtering process. 
Fortunately, such signals have a tendency to produce coherence maps with multiple peaks of 
similar strength (see the middle panel of Figure [Sll. We are able to reject this class of events by 
comparing the difference of the main peak of the coherence map to the second strongest peak. 
This quantity provides a cut value to reduce the probability of a weak CW signal passing as an 
event of interest. 

Anthropogenic noise is unpredictable and often involves weak CW signals in several bands. 
For this reason we optimize the level of the cut according to the amount of signal surviving 



it. Figure 12 shows the distribution of the difference between the first and second peak for 



calibration signals as well as simulated neutrino signals. The cut is set to preserve all the signals 
while removing a large portion of weak CW residuals. The thermal event distribution is also 
plotted showing a sharp exponential drop. We allow some thermal events to pass given that the 
linear discriminant, discussed in Section 3.2, since the linear discriminant already takes care of 
any potential thermal contamination in the data set. 
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Figure 9: Angular reconstruction error versus the single antenna impulse signal to noise ratio (SNR). The calibration 
pulser is a broadband source with most of its power at 200-400 MHz, falling off at higher frequencies. An exponential 
function is fitted to each trend. The angular reconstruction error in elevation (Theta) is generally better due to the fact 
that ANITA has long vertical baselines with antennas pointed in the same direction. The azimuthal reconstruction error 
(Phi) is worse due to the horizontal separation of the antennas being smaller (~ 1 meter) and the antennas not being all 
pointed in the same direction. 
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Figure 10: Plot of the peak coherence map value M(fmnv) versus the signal to noise ratio (SNR) of the coherently 
summed waveform VE(^f,„aA-). The distribution for 8 million simulated thermal noise events (blue contours and yellow 
dots) shows that, at the tail of the distribution, the peak coherence M(Xinax) and summed waveform SNR VY,(t, ?„„,-) 
are anti-correlated. A contour cut is selected that admits a thermal noise event with a probability of 0.01 for the whole 
ANITA-I flight. The red contours are the distribution of the whole ANITA-I flight data for the vertically polarized 
channels. The distribution is primarily thermal and matches the simulations well. The tail to the left (high M(Xmax) for 
low VY.it, T„,ax) SNR) is due to residual carrier wave signals. The contours that follow a diagonal (high M(f„,o.v) and high 
VY.(t,fmax) SNR ) are primarily due to anthropogenic signals. The black dots show triggered simulated neutrino events 
with energy of 3 X lO"* eV. These events are distributed along the diagonal of this plot as expected. 



The effect of weak CW interference also manifests itself in the comparison of the peak coher- 
ence value M(f„M v) versus the SNR of Vi,{t, r,nax) (the maximum peak-to-peak distance of the 
coherently summed waveform calculated in the direction of peak coherence). The population 
of events in the data that produce the contours on the region between an SNR of 1 and 3 and 
coherence values greater than 0.1, in Figure [TO] is due weak CW contamination, which results 
in an increased coherence while not displaying a strong SNR in V^C^ ^max)- The distribution of 



coherent waveform sum SNR is shown in Figure 13 for caUbration signals and simulated data. 
Cutting data with SNR<4 filters out most weak CW interference while retaining signals with 
100% efficiency. 
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Figure 1 1 : Histogram showing tlie linear discriminant from Figure |10| Tlie blue histogram represents thermal noise 
simulated data for 8M events, which is the number of events recorded with ANITA-I. The green line is an exponential fit 
to the tail of the distribution from simulations. The cut value is set at a value corresponding to 0.01 thermal background 
events leaking into the analysis. The red line shows the discriminant for the ANITA-I data for vertically polarized 
channels. 



4. Identification and Characterization of Weak Signals 

If a stationary source is too weak to be detected in a single interferometric image, the images 
can be stacked to provide a strong detection. This technique can be useful for characterizing 
known weak sources of noise as well as identifying new ones. We demonstrate the use of aver- 
aging interferometric images with short integration times to observe the Sun. We apply a similar 
technique to producing RF source maps of the Antarctic continent. 



4.1. Solar Imaging 

The Sun, while a significant thermal source in ANITA's frequency range, is not distinguishable 
in a single interferometric image created from ANITA data. Each ANITA interferometric image 
contains only ~100 ns of data. To detect such a source, one must combine many events by 
averaging the individual interferometric images together. Selection of a coordinate system in 
which the source of interest is stationary allows the signal to remain constant while the noise 
averages down, and is essential to this process. 

So far, we have plotted the coherence map M(f ) as a function of payload elevation Op and 
azimuth (pp. This is the coordinate system used in Figures l5] [6] [7] and [8] However, we can just 
as well select another reference coordinate system for the image. For the purpose of imaging the 
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Figure 12: Histogram showing the distribution for the difference between the first and second highest peaks of the 
coherence map M(f). The black histogram is for simulated thermal events and has a sharp exponential drop. The red 
histogram is for the calibration impulses at the Long Duration Balloon Facility in Williams Field, Antarctica. The green 
and blue histograms are for simulated neutrino signals with an energy of 3 X lO"* eV and lO" eV, respectively. The 
purpose of this cut is primarily meant to address carrier wave signals of anthropogenic origin, which tend to produce 
many peaks with comparable values. The magenta line allows a few thermal events to pass as signals while preserving 
signals of interest with 100% efficiency. When combined with the thermal linear discriminant cut the number of thermal 
events passing is negligible. 



Sun, we use solar azimuth 0s, where (ps - corresponds to the location of the Sun. Although 
at South Polar latitudes the Sun changes elevation angle throughout the day, it is reasonably 
stationary in the 30 minute time scales used for averaging. Making the coherence map in this 
coordinate system reduces to recalculating the baseline delays T(0p, (ps) and averaging the cross- 
correlation coefficients for each baseline Cij{T{6p, (ps)) to produce the coherence map M{6p, (ps)- 
Each map, using a ~100 ns snapshot in the case of ANITA, will not reveal an image of the Sun. 
However, the average map over A^ events 



1 ^ 



(13) 



where the index ev runs over all events used in the average, will produce an image given that A^ 
is large enough. 
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Figure 13: Histogram for tlie coherently summed waveform Ve(?, T,„ax) signal to noise ratio distribution of signals. The 
red and blue histograms are signals from the borehole and ground calibration pulsers, respectively, at the Long Duration 
Balloon facility in Williams Field, Antarctica. The green and red histograms are the distributions for simulated neutrino 
events at 3 X lO'* eV and lO" eV, respectively. The lack of events below an SNR of 4 is due to the ANITA trigger 
threshold. A cut on values of SNR < 4 filters most of the CW residual signals shown in FigurefTo] 



Figure 14 shows an image of the Sun S{6p, (ps) created by averaging together 10,000 ANITA 
events (~1 ms of data, recorded over a ~30 minutes period). The image created from the average 
amplitudes of each pixel clearly reveals the Sun and its surrounding sidelobes. This image shows 
a peak corresponding to the location of the Sun at Op ~20°. In addition, the image shows a bright 
spot consistent with the Sun's reflection on the ice at 9p — 20°. The reflection is stronger in the 
horizontally polarized image than the vertically polarized image, as expected from the Fresnel 
reflection coefficients on the surface of the ice, with index of refraction n - 1.35. Another 
interesting feature are the horizontal bright lines spanning the full azimuthal field of view. These 
lines could be due to the horizon where the noise temperature transitions from 270° K on the 
surface to 10° K on the sky. No conclusive evidence has been found to distinguish it from a 
detector artifact and the effect is currently under investigation. 

This technique is being further explored for a variety of applications and will be treated in 
detail in an upcoming publication fiE\. The Sun provides a constant source that allows us to 
monitor the antenna gain calibration throughout the flight. We can also monitor the surface radio 
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Figure 14: An image of the Sun formed by averaging 10,000 interferometric images of ANITA events. The location of 
the Sun is tracked throughout the three hours used to form the image and is set as a reference point. The Sun is at 20° in 
payload elevation angle. The azimuth angle coordinate is centered on the expected payload azimuthal location of the Sun. 
The fainter peak at -20° elevation is consistent with the reflection of the Sun on the ice. The reflection point is stronger in 
the horizontally polarized channels as expected. The color scales in the coherence map are normalized to the peak value 
of the map. Features immediately surrounding the image of the Sun are the sidelobes of ANlTA's point spread function. 
The straight bright curves, which extend across all azimuths, could be associated with the horizon (marked with a black 
line) as seen from ANITA's latitude of ~ 80° south at an altitude of 37 km. 



reflectivity by comparing the brightness of the Sun and its reflection both in the vertical and 
horizontal polarizations. This comparison provides a direct measure of the index of refraction of 
the surface of the ice. Surface roughness estimates are important to the energy determination of 
reflected ultra-high energy cosmic ray air shower events |26|. 



4.2. Man-Made RF Activity on the Antarctic Continent 

Except for a small fraction of signals expected from ultra-high energy particles, ANITA 
records two main types of events: first, thermal noise fluctuations that trigger the system compris- 
ing >99% of the data and second, anthropogenic signals originating from Antarctic bases, field 
camps, traverses, and potentially aircraft comprising the remaining <1% of the data. There is, 
however, some overlap in these events. Some anthropogenic sources can be well below ANITA's 
trigger threshold but still present in the thermal noise triggered data. Much like the treatment of 
the Sun in section 4.1, we are able to make subthreshold maps of the Antarctic continent. 
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For ANITA-I, a total of 8 M events were collected, making the projection of all coherence 
maps onto a coordinate system covering Antarctica impractical. However, we can take the peak 
value and direction of the coherence map M(f „,„ J for each event and project it onto the Antarctic 
continent. In the following, we describe the mapping procedure. First, we find the peak direction 
of the coherence map in payload coordinates 9p, ^p. We then transform the direction to a local 
East, North, Vertical (ENV) coordinate system to obtain the direction of the peak in Oenv, ^env- 
This involves correcting for payload heading and attitude offsets from vertical [TOl. The location 
of the payload is determined by an on-board global positioning system (GPS) unit reported in 
Earth-Centered Earth-Fixed (ECEF) Cartesian coordinates. The direction of the peak coherence 
also needs to be transformed to Oecef, 4>ecef- Given the position and direction of the peak 
coherence, we can propagate the ray to a location on the Antarctic continent. To do this, we bin 
the continent in Easting, Northing coordinateajin 10 km by 10 km squares. Using an elevation 
model of the continent |35 1 we determine the ECEF coordinates {xk,yk, Zk) for the center of each 
bin k. We then linearly propagate the ray from {xp,yp, Zp) in the direction of Oecef, 4'ecef to find 
the closest bin {xk,yk,Zk) consistent with that propagation. 



Figure 15 shows the average of the peak coherence values for events projected into each bin on 
the Antarctic continent for vertically polarized data of all 8 M events from ANITA-I. The blue 
background is consistent with the thermal noise expectation, while the colored regions indicate 



hot-spots of anthropogenic activity. Figure 16 shows a similar map made using all 21 .2 M events 
from ANITA-II |,29J . More anthropogenic hot-spots were observed with ANITA-II because of 
the increased exposure and sensitivity of that flight. 

When these maps, made for all the data, are compared to the event clusters passing the thermal 
noise and CW filters, we find clusters of activity that do not appear in the latter population of 
events. This indicates that there is potential for a single anthropogenic transient event that could 
be detected from among the background of sub-threshold signals. We therefore exclude any 
events from such sites in the final analysis |29|. It is also interesting to note that this technique 
provides a diagnostic for the radio quiet properties of a given site. If we look at how the events in 
each bin are distributed in peak coherence, we can determine whether the site is consistent with 
pure thermal noise or if there is some sub-threshold tail indicating weak anthropogenic activity. 
We have used this technique for visual inspection and we are currently looking into applying it 
more formally to the data analysis. 

Reliable directional reconstruction of anthropogenic activity was critical to the results of the 
ANITA-I and ANITA-II ultra-high energy neutrino searches If34ll29l . The criteria for tagging 
an above threshold event (an event passing the thermal noise and CW filters) as anthropogenic 
is any one of the following: the event clusters with other above threshold events, the event 
clusters with a known location of human activity regardless of whether it also clusters with 



other above threshold events, the event clusters with a local maximum (hot-spot) from Figure 15 



for ANITA-I or Figure 16 for ANITA-II which may be formed by purely subthreshold signals. 
If any of these cases are satisfied the event is rejected. For ANITA-I, one isolated event was 
identified as a man-made background event by association with a hot-spot formed solely out 
of subthreshold events. Even though there is a significant amount of human activity on the 
Antarctic continent, especially as seen by ANITA-II, only 36% of the continent was excluded 
in the ANITA-II neutrino search analysis due to proximity to anthropogenic events |29|. The 



^Easting and Northing coordinates place the South Pole as the origin with the y-axis (Northing) pointed along 0° 
longitude line and the x-axis (Easting) pointed along the 90° longitudinal line. 
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blue regions shown in Figures 15 and 16 indicate quiet regions where neutrino candidate signals 
could be found. 



5. Outlook and Conclusions 



The interferometric techniques developed for impulsive signals are applicable not only to event 
reconstruction but also to filtering of thermal noise and radio interference as well as the identi- 
fication of weak background sources. The methods presented are suited for antenna arrays with 
digital sampling capabilities, the kind of detector fit for searches of UHE particle impulsive tran- 
sient radio emission. The applications discussed in this paper are only a subset of the capabilities 
provided by the interferometric technique and several others are currently under development. 

Future applications of the Solar imaging technique include estimation of the surface roughness 
of the Antarctic ice by comparing the direct and reflected images. This technique can also be 
applied to a determination of the index of refraction and reflectivity of the Antarctic surface. The 
analysis of Solar reflections will also provide a measure of surface roughness on the Antarctic 
continent relevant to the energy determination of ultra-high energy cosmic ray events 1361 . 

The interferometric principles presented in this paper are also being applied towards the de- 
velopment of a new trigger for the third flight of ANITA. The interferometric trigger will contin- 
uously digitize data with 3-bit resolution in real-time. The beam-forming can then be performed 
in real-time using an application-specific integrated circuit or a field-programmable gate array, if 
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Figure 15: Coherence map of the Antarctic continent for the ANITA-I flight. The image is formed by obtaining the peak 
value of the coherence map of each ANITA event and then projecting its direction onto the Antarctic continent. The color 
scale indicates the average value of the peak of the coherence map M{r„,ax) for all triggered events which project onto a 
given bin on the map. This map has been useful for identifying the structure and location of anthropogenic backgrounds. 
Since no cuts are placed on the data prior to its projection onto the continent, it allows for the identification of regions of 
anthropogenic radiation that are below the threshold of the neutrino search cuts. Note that a large portion of the Antarctic 
continent (in dark blue) is radio quiet. McMurdo, on the bottom, and South Pole station, in the middle, are very strong 
emitters and dominate the radio anthropogenic background. 
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Figure 16: Coherence map of the Antarctic continent for the ANITA-II flight. The image is made in the same fashion 
as that shown in Figure [T3] McMurdo displays the largest and strongest emission. ANlTA-11 did not fly over the South 
Pole like ANITA-I but it was still able to observe strong anthropogenic backgrounds in its vicinity. The remaining bright 
spots are emissions from field camps and other radio experiments. 



the power consumption allows. The development of the interferometric technique for impulsive 
transients applied to hardware algorithms will be treated in a future publication. 

The application of interferometric imaging to impulsive signals has proved to be a powerful 
technique. It is by no means limited to what has been presented in this paper and it will continue 
to be exploited in future efforts for the radio detection of ultra-high energy particles. 
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Appendix A. Statistical Pointing Error Estimation 



Appendix A.l. Probability Density Function and Likelihood Estimator 

In this appendix we estimate the pointing errors from a likelihood estimator approach. Fol- 
lowing |33 1, the probability density function p(a, (p), for the amplitude a and phase (pof a phasor 
a = ae"^, resulting from the sum of a number of random contributions large enough to satisfy the 
central limit theorem, is 



p(a, (p)dad(p — 



Incr^ 



exp 



20-2 



dad<p. 



(A.l) 



where cr is the noise level. The probability density function has units of inverse amplitude times 
inverse radians and is written as p{a, (p)dad(p to make it explicit that it is a differential probability 
density. Expressed in terms of the real and imaginary parts of 5 = x + iy, the probability density 
function is a bivariate Gaussian distribution 



p(x, y)dxdy ■ 



1 



Incr^ 



exp 



x^ +y2 



2cr2 



dxdy. 



(A.2) 



In the presence of a signal phasor with s - s cos S + is sin 6 with amplitude s and phase S, in a 
noise background cr the probability density function becomes 



p(x,y)dxdy ■ 



1 



Incr^ 



exp 



(x- s cos <5)^ + (37 - ,? sin 6)'^ 



2cr^ 



dxdy. 



The likelihood function associated with the probability density function described is 

(x-scosdf + {y- ssind)^ 
- 2 log L(x, y; s, 6) = r 



(A.3) 



(A.4) 



Appendix A.2. Likelihood Estimator of a Digitized Signal Observed by an Array of Antennas 

Let's assume we have an array of Na antennas indexed by /. If the signal at each an- 
tenna is digitized with sampling frequency /, with Ns samples, then the frequency resolu- 
tion is A/ - fs/Ns- For uncorrelated noise, each frequency bin of width A/ is an inde- 
pendent measurement indexed by k. For an array of antennas with effective height /i,_t(0, 4>) 
observing an electric field Ei^ (in one polarization), the antenna signal can be represented as 
Si,k = Eithj^k(0, (f)) (cos a)iiT{6, (p) + ismajkT(9,(p)) where Ti{0,(p) is the geometrical delay of the 
signal cTj = R ■ f . For simplicity we have assumed that the overall electric field phase is at zero 
degrees and that the antenna effective height does not affect the phase other than by the geo- 
metric delay. Adding these dependencies is straightforward but we want to keep the following 
derivation as simple as possible. 



Using Equation A.4 the likelihood function for the array is given by 



N,/2 N^ 



-2 log L({xi,k}, {ya); E,co,0,cf>) ^^^Tj 



{xi^k - Ekhi^k cos cokTiiO, 4i)f + iyi^k - Euhj^k sin WiT,^^, ^)f 



k=0 i=\ 

(A.5) 
where {jc, _/.) and {y, j) are the measured real and imaginary parts for the phasors of each antenna 
indexed by / at each frequency indexed by k. The sum over frequencies has Ns/2+l independent 
contributions for a real digitized signal. 
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The error for is estimated from the second derivative of the maximum likelihood estimator. 
The first derivative with respect to Q gives 

d(-2 log L)/d0 = > > -^ — - [xi^k sin cotTi - y,;^ cos wtT,J (A.6) 

k=Q 1=1 i,k 

The second derivative with respect to 6 gives 



N,/2 Na 



dH-2iog lydo" = X Z ^^^-t 



k=0 1=1 "^i.i 



2 



Wi ( -^ (x,-,ii COS WiTi + yi^k sin wtT,) + -—^ (x,;^ sin w^t,- - y,- ^ cos wj^t,) 



(A.7) 
In the limit where the data approaches the modeled values x,^ — > Eiihjj. cos w^t, and yij^ — > 

^t/Za sin W^T; 

5V21ogL)/5^ =2f;|;(^^)' (A.8) 

Note that Ei^hij^lcri^k is the signal to noise ratio (snrij.) at the antenna / at frequency bin corre- 
sponding to COk. 

For the purpose of illustration, let us assume we have a collinear array so that t, = (Ri/c) cos 0,, 
where /?, is the distance to the origin and 0, is the angle between R and f . Let us also assume that 
only one frequency bin has stirik + 0. The expression then becomes 

d^(-2logL)/dd^ ^2j^isnn—Rsm6i\ (A.9) 

We estimate the angular error erg using the relation cTg'^ - {l/2)d^(-2logL/dd^). Note that the 
error is not arbitrarily small given the choice of a distant origin. If this is the case, all the angles 
0i will be small compensating for the choice of large /?,'s. Let us assume that a signal is incident 
orthogonal to the collinear array axis. If we set the origin at the location of one of the antennas, 
say the A^^th one, then its contribution drops out of the sum because R^a = 0. Then we have 

crs = , "^'" (A.10) 

Assuming that stiri = stir is the same for each antenna 

<rg^ '^1_ (A.11) 



•VSP^ 



For the ANITA geometry, the antenna beam pattern is wide enough so that there are several verti- 
cal baseline pairs contributing to the elevation error. In this case there is a fixed vertical distance 
giving an improvement of roughly ^INa - L The error can be calculated more rigorously using 
the effective heights of the antenna which make the channel to channel snr vary. 

In general, the digitization with A^, records is sensitive to A^s/2 H- 1 independent frequencies. 
Then we have 

A',,/2 Na , n x2 

d\-2\ogmdQ' = y y Lv«r,;i27r- sinfl; (A.12) 
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Let us assume that there is a relatively large number of frequency bins Nf, spanning indices k^ 
to ^1, for which the signal to noise ratio snrii. - snr is constant and non-zero. The error estimate 
is 



c 1 1 

<Te^-^ i ^= (A.13) 

Let us evaluate JYjk=k fk- ^^^ ^ digitized signal with A/ we have ft — kAf. For a large 
number of frequency bins the sum results in Yj^Lk ^^ ~ (^i ~ ^o)/3- Factorization gives k^—kQ = 
{ki — ko)(k^ + k^). The first term is the number of frequencies with non-zero SNR (ki - ko) = Nf. 
The second term gives Af^{k^ + k^ji^c^) ~ l/A^ for the central frequency A = 2cl({k\ + ko)Af) 
Together these give 

ere ^ ^-^ , ^ (A. 14) 



Thus, the direction error is reduced by a factor of yjNf. Intuitively, this result can be interpreted 
as each independent frequency providing an independent interferometric estimate of the incident 
angle 6. 

Appendix A.3. Likelihood Estimator of an Interferometric Array 

A similar likelihood analysis can be performed on the cross-correlation phasor defined as the 
product a,a* of two phasors. The cross-correlation phasor for a signal incident on the inter- 
ferometric array is given by E^hjkhjk cos w,t(Tj - tj) + i sin cokiri - tM. The advantage of this 
approach is that any ambiguities due to choice of coordinate system vanish, since only delay 
diff'erences between antenna pairs are counted. However, the estimation of the cross-correlation 
phasor noise is non-trivial and the accounting of independent contributions has additional com- 
plications. A maximum likelihood analysis technique along these lines has been developed for 
interferometric observations of the cosmic microwave background 1371 . This approach is cur- 
rently being developed as a potential improvement on the interferometric analysis presented on 
this paper and for a more rigorous accounting of errors. This will be the subject of a future 
pubhcation. 
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